Noise Correlations in Three- Terminal Diffusive Superconductor-Normal 
Metal-Superconductor Nanostructures 
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We present measurements of current noise and cross-correlations in three-terminal 
Superconductor-Normal metal-Superconductor (S-N-S) nanostructures that are potential solid-state 
entanglers thanks to Andreev reflections at the N-S interfaces. The noise correlation measurements 
spanned from the regime where electron-electron interactions are relevant to the regime of Inco- 
herent Multiple Andreev Reflection (IMAR). In the latter regime, negative cross-correlations are 
observed in samples with closely-spaced junctions. 
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Non-local entanglement - the emblematic ingredient to 
quantum physics - has been proposed and debated for 
nearly more than a century [jy]. Having already been 
demonstrated in photons [2], it remains difficult to ob- 
serve for massive particles like electrons. One interest- 
ing approach is to perform noise-correlation measure- 
ments using superconductors, as they are natural sources 
of entangled electrons. In mesoscopic systems, non- 
equilibrium current noise measurements provide informa- 
tion on the charge and statistics of current-carrying states 
0. For non-interacting electrons, the Pauli exclusion 
principle dictates the zero- frequency cross-correlations to 
be negative Q , whereas no such rule applies in the pres- 
ence of interactions. In early experiments on electronic 
analogues of teh Hanbury-Brown and Twiss jSj experi- 
ment, noise correlations came out to be negative. For 
beams of partitioned electrons using high mobility GaAs 
two-dimensional electron gases (2DEG), anti-correlations 
were revealed, thus illustrating that fermions exclude 
each other . Fermionic correlations have also been ob- 
served for free electrons and neutrons Q. Interest- 
ingly, positive correlations observed in a purely normal 
2DEG were ascribed to different scattering mechanisms 
within the device 9;]. Recently, cross-correlations with 
a bias-dependent sign were reported in a three-terminal 
superconducting hybrid nanostructure with tunnel con- 
tacts p^]. 

In a superconducting hybrid (N-S) beam splitter made 
of two normal metal (N) leads in contact with a super- 
conductor (S), an incident electron (hole) from one lead 
can be reflected at the superconductor interface as a hole 
(electron) propagating into the other lead. This non-local 
crossed Andreev reflection (CAR) process corresponds to 
the transfer (creation) of a superconducting Cooper pai r 
into (from) two entangled electrons in the two leads [Uj . 
It creates positive correlations between currents flowing 
in each of the two leads [12] . On the contrary, elastic 
co-tunneling of an electron (hole) from one lead to the 
other contributes to negative correlations [l^l- Whereas 
negative correlations are expected from a semi-classical 



approach neglecting proximity effect [l4| , positive corre- 
lations at sub-gap energies are predicted for intermediate 
values of transparencies [15]. Positive correlations could 
also arise from synchronized Andreev reflections [l6l |. In 
a normal metallic dot connected to all superconducting 
leads, positive or negative correlations are expected, de- 
pending ripon properties of contacts between the dot and 



its leads 17 



In a S-N-S junction, charge transfer is mediated by 
Multiple Andreev Reflections (MAR), a process during 
which quasiparticles undergo successive Andreev reflec- 
tions at both interfaces until their energy reaches the su- 
perconducting gap. At a bias above the Thouless energy 
Exh ~ TiD/L^, where D is the diffusion constant and 
L is the junction length, MAR are not phase-correlated. 
In this incoherent MAR (IMAR) regime, the noise has 
been found to be very much enhanced compared to the 
normal case due to the confinement of the sub-gap elec- 
trons in the sandwiched normal metal 18|. In the same 
regime, cross-correlations are also predicted to be en- 
hanced due to IMAR processes [17|. At low bias voltage. 



when the time taken by the quasiparticles to reach the 
superconducting gap exceeds the electron-electron inter- 
action time, the IMAR cycle is interrupted. In this hot 
electron regime, a Fermi-Dirac like distribution with an 
elevated effective temperature is then restored ^1^ . 

In this Letter, we present measurements of noise 
and cross-correlations in three-terminal diffusive S-N-S 
nanostructures. In the IMAR regime, we demonstrate 
negative correlations in samples where the two junctions 
are closely spaced. 

Samples comprising three-terminal double S-N-S junc- 
tions (see SEM pictures in Fig. [T]) were fabricated by mul- 
tiple angle evaporation through a PMMA-PMMA/MAA 
bilayer mask in an ultra-high vacuum chamber. Evap- 
oration of Cu with a 50 nm thickness was followed im- 
mediately by the evaporation of 500 nm-thick Al elec- 
trodes, thereby forming diffusive S-N-S junctions across 
the Cu bridge with highly transparent S-N interfaces. 
The width of the Cu part was 0.9 /im. The samples have 
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FIG. 1. Current-voltage characteristics of each junction Ji 
and J2 of samples W and C at 100 mK. Insets: SEM pictures 
of both sample geometries. The distance between junctions 
is 10 fim for sample W and 0.5 fim for sample C. 

two different distances between junctions: 10 /im for the 
widely-separated junctions (sample W) and 0.5 for 
the closely-spaced ones (sample C). In the latter case, 
their distance is comparable with the superconducting 
Al coherence length, thus making CAR and elastic co- 
tunneling probable. Here, we discuss results obtained on 
two samples W and C with a junction length L = 1.3 fim. 

Figure [T] depicts I-V characteristics of individual junc- 
tions of the two samples at 100 mK. The junction re- 
sistances of a given sample are very much symmetrical, 
although they differ slightly in the widely-separated ge- 
ometry. All junctions show a superconducting branch 
with a small critical current (of the order of few /xA) , fol- 
lowed by a linear part corresponding to the normal-state 
resistance of the normal metal, of the order of 1.5 J7. The 
related Thouless energy is estimated to be about 5 /ieV. 
The abrupt transition in I-V curves at large currents (> 
50 /iA) is due to depairing effects in the superconducting 
Al electrodes. The temperature dependence of the I-V 
characteristics (not shown) does not show much of vari- 
ation, apart from the depairing current, which decreases 
at higher temperatures and vanishing critical currents. 

We have used a new experimental set-up especially de- 
signed to measure current fluctuations and noise cross- 
correlations in three-terminal devices at low tempera- 
tures. The experiment operates down to 30 mK and 
is equipped with three commercial SQUIDs (Supercon- 
ducting Quantum Interference Devices) as sensitive cur- 
rent amplifiers, see Fig. [2] Each junction (arm) of the 
sample is connected to the input coil of a SQUID sitting 
in the Helium bath at 4.2 K. A reference resistor Rref 
of low resistance (0.092 il) together with a third SQUID 
coil is connected in parallel to the sample for voltage bi- 
asing. Voltage probes allow us to measure the voltage 




FIG. 2. Schematics of the circuit consisting of the sample, 
each junction being represented by differential resistors Ri 
and i?2, and the three SQUIDS to measure the current fluc- 
tuations. The reference resistor Rref is used to voltage-bias 
the sample. According to Nyquist representation, each re- 
sistance in the model is associated with a current source in 
parallel. The two additional resistors r' and r" (r = r' + r") 
are due to contact between distinct superconducting elements. 
Their temperatures are the bath temperature and 4.2 K re- 
spectively. Their associated noise source is not shown for 
clarity. 

drop across each junction (Vi, V2) and across the refer- 
ence resistor (Ke/). We found that Vi ~ V2 but Ke/ 
differs slightly from Ui(2). This difference, of about 5 
%, is due to an additional resistor between the sample 
and the reference resistor. For the two samples discussed 
here, it corresponds to a resistance r ~ 40 mf2, which is 
not negligible compared to Rref and needs to be taken 
into account in the model discussed below. The intrinsic 
noise level of each SQUID expressed in equivalent current 
at its input coil is of few pA/-\/Hz. At frequencies above 
few hundreds of Hertz, 1// noise contributions were neg- 
ligible for all bias currents, so that a frequency domain 
between 800 Hz and 4 kHz was chosen for the noise and 
correlation measurements. 

In accordance with Nyquist representation, each resis- 
tive part of the circuit is associated with a current source 
Sli. The two junctions of a sample are represented as 
two non- linear resistors i?i and i?2. Each SQUID mea- 
sures partially each of the 3 current sources [19|. From 
the 3 fluctuating SQUID currents SIl^, we can perform 3 
auto-correlations ACi = 511^511^ and 3 cross-correlations 
XCij = 511^611^ non-independent measurements. The 
spectral densities of noise ACi and correlations XCij 
are related to the physical quantities Si = SliSIi and 
Sij = SliSIj through a 6 x 6 matrix. Here, and iS'2 are 
the noise of each junction of a sample and S3 the thermal 
noise of the resistor {Rref + r). 

Figure [3] shows the 6 raw datas {ACi} and {XCij} as 
a function of the voltage drop across sample W. The sum 
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FIG. 3. Raw data of Ad's and Xdj's from sample W. The 
sum represented by the curve along about 20 pA^ /Hz is con- 
stant, as expected, but non-zero due to the set-up background 
noise. 



ACi + AC2 + AC3 + 2{XCi2 + XC13 + XC23) is constant 
thanks to current conservation law, but non-zero due to 
the noise and correlation backgrounds of the experimen- 
tal set-up. In order to extract the quantities of interest 
Si , iS'2 and 5*12 , one needs to choose 3 independent mea- 
surements. We chose to focus on ACi, AC2 and XC12, 
which, up to the first order, are close to Si, S2 and S'12. 
In a matrix form, the equations of the system reduce to: 



512 



■ ACi - AC° ' 

AC2 - AC^ 
XC12 — XC12 



[N] S3 (1) 



where ACi, AC2 and XC12 represent the background 
noise and correlation of the set-up. Finding Si, S2 and 
S'12 thus relies on the determination of the noise and cor- 
relation backgrounds, as well as of the thermal noise S3. 
The elements of the 3x3 matrix [M~^] and of the vec- 
tor [TV] depend only on the differential resistance of each 
junction and on the resistance Rref + ^■ 

Let us first discuss the case of sample W, where zero 
correlations (S'12 = 0) are expected. For a S-N-S junction 
in the IMAR regime, the shot noise in the diffusive case 
2e//3 is amplified by the number 1 -I- 2A/eV of times a 
quasi-particle is Andreev-refiected before it reaches the 
gap edge [23|: 



1 2A 

SiV)^-2eI[l + -] 



(2) 



where A is the superconduc ting gap. Neglecting the con- 
tribution of XC12 to Si(2) [l9|, we adjusted ACi, AC2 
and A = 170 ^eV so that Si and S2 fit this expres- 
sion. As inelastic collisions are neglected here, Eq. [5] is 
valid in a voltage range that is limited but enough to 



FIG. 4. Spectral density of noise Si and 5*2 and cross- 
correlation 512 of sample W at 100 mK, obtained using Eq. 
[1] The lines correspond to the IMAR noise prediction Eq. 
|2] for each individual junction. The dashed region indicates 
the interacting hot electron regime, where Eq. [2] is no longer 
valid. 



achieve a reliable fit. Knowing AC^ and AC2 and choos- 
ing XC12 — first, we can plot the correlation S12. Our 
analysis shows that, in order to obtain a bias-independent 
S12, we need to consider the spurious resistance r as be- 
ing split into two parts r' and r" sitting at the mixing 
chamber with a temperature T and 4.2 K respectively. 
The best result is obtained for r' = 15 mfi and r" = 25 
mfl, see Fig. |4l 

Finally, by choosing XC12, we get S12 equal to zero 
over almost the entire voltage range. It differs from zero 
only at higher voltages approaching the dcpairing regime, 
where the common superconducting electrode becomes 
resistive and the model is no longer valid. From cali- 
bration measurements, we know that, due to the current 
bias, the electronic temperature of Rref can reach up to 
300 mK from a bath temperature of 100 mK. However, 
this has only a minor effect here since most of the noise 
S3 comes from the resistance r" sitting at 4.2 K [l9| . 

For sample C with closely-spaced junctions, the mea- 
sured contact resistance r has the same value as that 
of sample W. The results can therefore be analyzed ex- 
actly the same way as above. In order to adjust the 
correlation background XC12 value, we consider the low 
bias regime. Here, the noise of each junction is thermal 
with an elevated effective electron temperature. The two 
noise sources 5/i(2) can then be assumed as uncorrelated 
and hence the cross-correlation noise is zero at low bias: 
Si2(V^ ~ 0) = 0. We used this criteria to determine 
XC°2. 

The overall results of noise and correlations for the 
two different geometries are depicted in Fig. [5] for dif- 
ferent temperatures from 100 to 700 mK. Clearly, the 
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FIG. 5. Spectral density of noise Si and cross-correlation S12 
of the two samples W and C for various temperatures from 
100 to 700 mK. The noise 52 would give very similar results to 
that of S\. The solid line reproduces the IMAR predictions 
of Eq. [2] whereas the dashed line enlightens that negative 
correlations appear in the IMAR regime. 



two samples show different cross-correlation behaviors 
above a voltage of about 40 /iV, which corresponds to the 
cross-over between the hot electron regime and the IMAR 
regime. For sample C in the IMAR regime, correlations 
are negative, up to the depairing regime. When the tem- 
perature is increased, thermally activated quasiparticles 
generate additional thermal noise ^18., 2^ ■ As these fluc- 
tuations are uncorrelated, correlations are expected to 
exhibit a negligible temperature dependence. This con- 
trast in behavior is actually observed in Fig. [S] In ad- 
dition to the two samples discussed here, other samples 
with the same geometry but a different junction length 
of 1 fim were studied and produced very similar results. 
For these shorter junctions, the noise level is larger since 
the resistance is lower, see Eq. [2j Again, correlations 
were found to be zero in the sample with widely sepa- 
rated junctions and negative in that with closely spaced 
junctions. 

Negative correlations are expected for fermionic sys- 
tems. Therefore, our results suggest that negative cor- 
relations arise from partition of quasiparticles injected 
above the superconducting gap as the result of IMAR 
processes. When IMAR processes are interrupted by in- 
elastic collisions, the quasiparticle current is reduced and 
correlations vanish. This is in agreement with our find- 
ings. A quasiparticle current is known to vanish over the 
quasiparticle diffusion length, of the order of few fj,m in 
Al 



restrict the possibility to observe positive correlations in 
three-terminal devices with all superconducting contacts. 

In conclusion, we have measured negative cross- 
correlations in three-terminal diffusive S-N-S nanostruc- 
tures in the incoherent multiple Andreev reflections 
regime. This experiment opens the way towards a better 
understanding of non-locality and entanglement in su- 
perconducting nano-devices with various interface trans- 
parencies and bias schemes. This work was supported 
by ANR ELEC-EPR project. Samples were fabricated 
using NanoFab facility at Neel Institute. We thank M. 
Houzet, J. Meyer and R. Melin for discussions. 



21|. Zero correlations are thus expected in samples 
with widely separated junctions, as observed. To our 
knowledge, the role of the quasiparticle current on corre- 
lations has never been investigated theoretically. It may 
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Supplemental Material. 
SQUID currents equations. 

In accordance with Nyquist representation, each resis- 
tive part of the circuit is associated with a current source 
Sli with i = 1,2 or 3. The two junctions of a sample 
are represented as two non-hnear resistors Ri and R2, 
see Fig. 2 of the main text. The current fluctuations 
through each SQUID can thus be expressed as: 



611 



Ri 



Ri + Rjk 



R 



Ri + Rjk 



iSI,+6h) (3) 



whore 2, j and k G {1, 2, 3} and Rjk is the equivalent re- 
sistance of Rj in parallel to Rk ■ The third resistor R3 in 
the model correspond to Rref + r. From the three fluctu- 
ating SQUID currents <5J|q, we can perform three auto- 
correlations ACi = SllqSIlq and three cross-correlations 
XCij = dll^dljg non-independent measurements given 
by: 



ACi =< FFT*{dlig)FFT{6Ilg) > 
XCij =< FFT*{6I%)FFT{5li^) > 



(4) 



where FFT stands for the FastFourierTransform, 
FFT* its complex conjugate and < ... > the rms av- 
erage at the spectrum analyzer. 

The spectral densities of noise {ACi} and cross- 
correlations {XCij} are related to the physical quantities 
Si = SliSIi and Sij = SliSIj through a 6 x 6 matrix: 



' ACi ' 






AC2 




S2 


AC3 
XC12 




S3 
S12 


XC23 




S23 






Sl3 



(5) 



However, since Rref and r are macroscopic resistors sit- 
ting far apart from the sample, we consider that S'13 = 
S23 = 0. The SQUID current equations reduce then 
to a 6x4 matrix. As an example, the 6x4 matrix for 
R^=R^ = 1.5f7 and R3 = Rref +r = O.IQ is : 



[M6X4] 



0.885813 
0.00346021 

0.778547 
-0.0553633 
0.0519031 
, -0.83045 



0.00346021 
0.885813 
0.778547 

-0.0553633 
-0.83045 

0.0519031 



0.00346021 
0.00346021 
0.0138408 
0.00346021 
-0.00692042 
-0.00692042 



-0.110727- 
-0.110727 
1.55709 
0.889273 
-0.778547 
-0.778547. 

(6) 



As the six equations relating the six possible measure- 
ments, {ACi} and {XCij}, to the four quantities, {Si} 
and 512, are not independent, the set of equations needs 
to be reduced to three independent equations. For this, 
we assume that the spectral density of noise S3 gener- 
ated by the resistance R3 = Rref -|- r is only thermal. 
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The SQUID equations reduce then to : 



[N] S3 



' ACi' 




'Si' 




= [M] 


S2 


_XCi2_ 




_Sl2_ 



■ ACf ' 

ACl 
XCI2 



(7) 



where we have introduced the noise and correlation back- 
grounds AC?, AC^ and XC^^ of SQUIDS 1 and 2. The 
matrix [M] is now a 3x3 matrix and [N] a vector, whose 
elements depend only on the differential resistance of 
each junction and on the resistance Rref + r. 

To obtain the noise spectral densities Si and ^2 and the 
correlations S'12 we need to inverse the new 3x3 matrix. 
We then obtain : 



Si 
S2 
S12 



■ ACi - AC? ■ 

AC2 - AC^ 
XC12 — XC12 



\N\ Ss 



Again, as an example, the inverse matrix [MY 
Ri ~ R2 — 1.5J7 and i?3 = Rref + r = {)AVt reads : 



1.13778 0.00444444 0.142222 
0.00444444 1.13778 0.142222 
0.0711111 0.0711111 1.14222 



(8) 



for 



(9) 



Effect of the temperature of the spurious 
resistor r 

As the additional resistor r is due to contact between 
distinct superconducting wiring elements, it can be split 
into two parts r' and r" sitting at the mixing chamber 
with a temperature T and 4.2 K respectively. The noise 
power ^3 is then 

53 = 4fcs[(i?re/ +r)T + r" x 4.2]/(i?,e/ + rf (10) 

To analyze the noise and correlations results of sample 
W, we proceed in two steps. As a first step, we only 
consider 5*1 and ^2 and suppose that Rref and r are at 
the bath temperature T (r^ — AOmn and r" = in Eq 
Uni). We adjust ACf, AC^ and A = 170^eF so that 
Si and 52 fit the IMAR prediction (Eq. 2 of the main 
text), see Fig. 4 of the main text. Here, we neglect the 



contribution of XC12 to 5i(2), which is supported by the 
facts that XC12 is much smaller than ACi^2) and that 
the matrix elements coupling XC12 to 5i(2) are small 
numbers. 

Once we have estimated AC?, AC2, we can plot 5i2 
using equation [5] and choosing XC12 = first. At this 
stage, we see in Fig. |6] of this Supplemental Material, 
that the noise correlations 5i2 is not constant as a func- 
tion of the voltage (up triangle black symbols). We can 
also consider that the resistor r is entirely due to contact 
at the SQUID having therefore a temperature of 4.2K 
[r' = and r" = 40mf2 in Eq[TU]). Again the 5i2 is not 
constant (down triangle green symbols). We then adjust 



N 



1 5 m n @ 0.1 K and 25 mfl @ 4.2 K 
XC12 A 40mri@0.1K 

40 m £1 @ 4.2K 




20 40 60 80 
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FIG. 6. Cross-correlation S12 between the two junctions of 
the sample W at 100 mK for various choices of the additional 
contact resistance split between the mixing chamber and the 
4.2 K stages on the correlation amplitude Si2- 



r' and r" so that the noise correlations 5i2 becomes bias- 
independent. The new expression for 53 is then fed back 
in Eq. [S]for iterative and fine adjustment of the param- 
eters values. The best result is obtained for r' = 15 mfi 
and r" = 25 mfi. Finally, XC12 is adjusted to get 5i2 
equal to zero almost over the entire voltage range (square 
red symbols) . Since the contact resistance is found to be 
the same for sample C, we use the same split of the ad- 
ditional resistor r to analyze the data obtained on that 
sample with closely-spaced junctions. 



